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We study the coherent dynamics of one- and two-electron transport in a linear array of tunnel- 
coupled quantum dots. We find that this system exhibits a rich variety of coherent phenomena, 
ranging from electron wavepacket propagation and interference to two-particle bonding and en- 
tanglement. Our studies, apart from their relevance to the exploration of quantum dynamics and 
transport in periodic structures, are also aimed at possible applications in future quantum compu- 
tation schemes. 
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I. INTRODUCTION 

Quantum transport is by now an almost generic term 
encompassing a broad class of phenomena in normally 
quite distinct areas of physics. Traditionally the term was 
found in solid state physics 0, while through develop- 
ments over the last several years, it is commonly found in 
the context of nonlinear dynamics pertaining to quantum 
chaos , in cold atoms in optical lattices U , Rydberg 
atoms strongly driven by microwave fields or even 
laser-driven multilevel systems _6, _7., ,8J, the latter dat- 
ing back twenty five years or so, although the term was 
never used in that context. Issues of coherent propagation 
of wavepackets, their interplay with disorder and/or fluc- 
tuations, localization, etc. are now discussed in many of 
these contexts and often interchangeably f§|. 

A new context, namely low-dimensional nanostruc- 
tures and in particular quantum dots (QD) |^ ^ 
has emerged over the last ten or so years which ap- 
pears to overlap with and in some sense unify sev- 
eral of these areas. Often referred to as artificial atoms 
[l^ Il3| , they offer the unprecedented possibility to con- 
struct at will and explore situations ranging from prac- 
tically single atom to a fully solid state many-body sys- 
tems Hi El El 111 E El S El lii ^ ^j. The 

nanofabrication possibilities of tailoring structures to de- 
sired geometries and specifications, and controlling the 
number and mobility of electrons confined within a re- 
gion of space 0, 13, [2I , are some of the features 
that make these structures unique tools for the study of 
a variety of preselected set of phenomena, not to mention 
the potential for devices. Given the controllable quantum 
properties of the electrons in such confined structures, 
the possibility of their application to schemes of quantum 
computers [27ll23 | has not escaped attention as attested 
by an ev er growing ple thora of papers on this general 

issue IM Ells ImIm 111 111 lai. 
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Motivated in part by the connection with quantum 
computing and in part by the desire to explore the dy- 
namics of wavepackets in periodic structures, we began 
with the study of such wavepackets in the simplest case 
of one electron tunneling through an array of QDs. As 
discussed in some detail below, the problem turns out to 
have unexpected similarities with an old problem, namely 
the molecular excitation through a ladder of vibrational 
levels .Models developed in that context more 

than twenty five years ago found no application, sim- 
ply because a vibrational structure is too complex to be 
amenable to simple models of, for example, equidistant 
levels or exact levels of a harmonic oscillator. On the 
other hand, in the present case, one can easily contem- 
plate fabricating almost any desired level and tunneling 
rates arrangement. Thus, models that were of only aca- 
demic interest then, now become of practical relevance. 
At the same time, the equations governing the propaga- 
tion in a chain of quantum dots are similar to those found 
in the 1-dimensional (ID) Ising model js^. 

Even a slight generalization, namely the presence of 
one more excess electron in the chain presents an entirely 
different situation. Its counterpart in multilevel systems 
does not exist, while it is now the 2D Ising model (111 that 
is related to the chain of QDs. In our case, however, it is 
the dynamical rather than the thermodynamic properties 
of the system that we explore. With two electrons avail- 
able for tunneling, an additional set of phenomena such 
as Coulomb blockade, bonding (to be described below), 
and correlation come into play. Surprisingly, we find 
that preservation of coherence found in the one-electron 
case, under a judiciously chosen set of parameters, car- 
ries through in the two-electron case. It is this finding 
that brings in the possible connection with schemes for 
the quantum computer. Such a possibility could be con- 
templated when the electron spin is brought into the pic- 
ture, combined with coherent propagation within a chain 
of QDs and manipulation of spin-entanglement (29i i3Qj | , 
as discussed in Sec. [V] 

Even with the high degree of accuracy in nanofabri- 
cation technology, issues of randomness and fiuctuations 
0,^, originating from more than one physical effect, 
are an inevitable part of the picture. Their presence, in ad- 
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dition to requiring attention with respect to their impact 
on coherent propagation and entanglement, represents a 
class of further phenomena related to the well-known lo- 
calization induced by disorder, and related issues 0.We 
have therefore examined some of our basic results in the 
presence of randomness. Although not necessary for the 
main arguments of this paper, for the sake of complete- 
ness, we have also explicitly included in our analysis the 
effect of measurement. This has necessitated the adop- 
tion of Monte Carlo techniques which are well established 
in quantum optics and provide a direct account of the ef- 
fect of a measurement on the system [s^, IS '39] . 

We would like to state at the outset that, although we 
do realize that each quantum dot in general has many 
electrons, it is nevertheless possible to arrange the condi- 
tions so that only one electron per dot tunnels, while the 
rest remain "frozen" .It ma y n ot be routine, but on the 
basis of existing literature [III ITsL IT^ ^| , it is techno- 
logically possible. In fact many papers with proposals for 
quantu m reg isters employing QDs are based on that very 
idea [2^ |23, 113, ■ We do also realize that assuming a 
chain of 10 or 20 dots, as we do in the sections to follow, 
may indeed stretch present day technological capabilities. 
We do so, however, in the interest of exploring a number 
of unusual effects which if not quite feasible today will 
most likely be feasible tomorrow. But when it comes to 
entanglement in connection with a quantum register, in 
the last section of the paper, even two dots fulfill our 
requirements. 

We have structured this paper as follows: In Sec. ^ 
we present the mathematical formalism describing ID 
chains of QDs, upon which we build the theory of co- 
herent single-electron propagation in Sec. IIIII and two- 
electron propagation in Sec. IIVI In Sec. 0we show a 
scheme for entangling two electrons in the chain via their 
controlled collision. Our conclusions are summarized in 
Sec. ED 

II. MATHEMATICAL FORMALISM 

The system under consideration consists of a linear ar- 
ray of A'' nearly identical QDs which are electrostatically 
defined in a two-dimensional electron gas by means of 
metallic gates on top of a semiconductor heterostructure, 
e.g., GaAs/AlGaAs |l3,|Il[l3- We describe the evolu- 
tion of the system using the second-quantized, extended 
Mott-Hubbard Hamiltonian [T7Ll20ll2^ 

^ = H ""^ ^ ^ Unjirij - 1) 

(aLoja + H.c.) + ^ V^jn.rij, (1) 

i<j,a i<j 

the first two terms being responsible for the intradot ef- 
fects and the last two terms describing the interdot in- 
teractions. Here al^ and the creation and an- 
nihilation operators for an electron in state a with the 




FIG. 1: Schematic drawing of the chain of tunnel-coupled 
QDs, and the last dot is dissipatively coupled to an external 
single electron detector (SED). 



single-particle energy Sja and electronic orbital V'ia(r), 

U = 8i^/drdr'|V',a(r)nV.,„.(r')|V|r-r'| eV^,, 
with Cg ~ 8er£oR being the self-capacitance for 2D disk- 
shaped QD {er — 13 for GaAs), is the on-site Coulomb 
repulsion, rij = '^la'^ja is the total electron number 

operator of the jth dot, ty,^ = ^ J dr'ip*^{r)V'^'ipja{r), 
with m* being the electron effective mass (m* ~ 0.067me 
in GaAs) , are the coherent tunnel matrix elements which 
are given by the overlap of the electronic orbitals of ad- 
jacent dots {j = i + 1) and thereby can be controlled 
by the external voltage applied to the gates defining the 
corresponding tunneling barriers between the dots, and 
V^j ~ U{CICg)\'-^\, with C < Cg being the interdot ca- 
pacitance, describe the interdot electrostatic interaction. 
Note that in general, the index a refers to both orbital 
and spin state of an electron. In the Coulomb blockade 
and tight-binding regime, when the on-site Coulomb re- 
pulsion and single-particle level-spacing Ae ~ ^''r? 
much larger than the tunneling rates, U > /S.e ^ tij^a 
[4ll |. only the equivalent states of the neighboring dots 
are tunnel- coupled to each other. Therefore throughout 
this paper we will consider a single doubly- (spin-) de- 
generate level per dot (a G {Ti ■!}); assuming further that 
the tunneling rates do not depend on the electron spin 
itij,a = Uj)- In our subsequent discussion, we assume 
that the matrix elements Vij = V are non-vanishing for 
nearest neighbors only, j = i + 1. Interdot repulsion 
can be further suppressed in the presence of a nearby 
conducting backgate, where image charges of excess elec- 
trons completely screen the interdot Coulomb repulsion, 
in which case V ^ 0. 

We study situations where one or more preselected 
QDs are initially doped with single electrons (Fig. 
The preparation of such initial conditions can be ac- 
complished by first applying a large negative voltage to 
the gates defining the electron confining potentials, while 
keeping the tunneling barriers low, thereby depleting the 
chain, and then doping one or more preselected QDs with 
single electrons by lowering the confining potentials and 
carefully manipulating the tunneling barriers between 
the dots and the external electron reservoir |2^ |2^ . Fi- 
nally, the system can be isolated by closing the tunnelings 
between the two ends of the chain and the reservoir. To 
start the evolution, one then tunes the interdot tunneling 
rates to the preselected values. This process should, on 
the one hand, be fast enough on the time scale of t~-^ 
so that no appreciable change in the initial state of the 
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system occurs during the switching time Tsw, and, on the 
other hand, adiabatic so that no nonresonant coupling 
between the dots is induced: e~^, U^^ < Tgw < t'^j' ■ 

To examine the influence of disorder due to structure 
imperfections, gate voltage fluctuations [TH IT^ li^. 
electron-phonon interactions [s^ . etc., we allow for ran- 
dom fluctuations of the intradot energy levels and inter- 
dot couplings, which would result in a decoherence of the 
electron wavepacket propagation. Specifically, we model 
the energy levels Sj and the interdot couplings tij as 
Gaussian random variables, with mean values Eq, Iq and 
variances de,dt < Iq, respectively. To account also for a 
possible measurement process, we assume the last dot of 
the chain to be dissipatively coupled to an external sin- 
gle electron detector (SED) with a rate 7^ = 7 (Fig. 
Such detection can be realized by coupling the last dot to 
another larger measuring dot at lower potential, which is 
continuously monitored by a nearby quantum point con- 
tact conductor, realizin g a broad-bandwidth single elec- 
tron transistor puJi lisL ISj. Then, the electron tunnels 
out of the chain into an effective continuum represented 
by many available, broad (due to the continuous moni- 
toring), empty levels of the measuring dot. 

The treatment of this type of problems is amenable to 
the density operator approach [23, 113 - In the case of 
only one excess electron present in the system, one has 
to solve iV^ coupled differential equations (equivalent to 
the optical Bloch equations) for the density matrix ele- 
ments. When, however, the initial number Ue of excess 
electrons in the system is larger than one, N ^ > 1, 
the number of coupled density matrix equations to be 
solved scales approximately as iV^"=, which is compu- 
tationally demanding. Having in mind the subsequent 
application of the method to a chain of QDs doped ini- 
tially with more than one electron (see below), we adopt 
here an alternative but equivalent approach based on the 
Monte Carlo stochastic wavefunctions, which has been 
well developed and widely employed in the context of 
quantum optics Is^jlH,!!^. Thus, for a set of Gaussian 
random numbers Sj, Uj we propagate the wavefunction 
of the system |5'(t)) using the effective non-hermitian 
Hamiltonian Hgfi — H — ^"fnpj, where H is given by 
Eq. This involves solving ~ N"^ amplitude equa- 

tions. The fact that H^g is non-hermitian implies that 
the norm of the wavefunction is not conserved. The prop- 
agation is interrupted by quantum jumps (corresponding 
to detector clicks in a real experiment) that occur when 
||^(r)|p = r, where r is a uniformly distributed random 
number, between and 1. Thus, a quantum jump corre- 
sponds to the detection of one electron by the detector, 
while the loss of the electron from the system projects the 
initial subspace of N QDs sharing Ue electrons, to a sub- 
space of N QDs sharing rie — 1 electrons. The post-jump 
wave function of the system is thus given by 



" \/(*|aLaJVa|*) 
where the denominator ensures its renormalization. A 



new random number r is then generated which deter- 
mines the subsequent jump event in our simulations. The 
propagation is terminated once the system has reached 
the zero electron state |0) = |0i, ...Oat). This procedure 
is repeated and the results are averaged over a large num- 
ber of independent realizations (trajectories). Thus this 
method allows one to model both a single experiment 
outcome, corresponding to a single trajectory, as well 
as the ensemble average measurement, corresponding to 
the averaging over many independent trajectories, while 
the density matrix approach is capable to adequately de- 
scribe only the latter. Since in the quantum computation 
schemes, one usually deals with a single measurement on 
the qubit [23| (with the notable exception of NMR-based 
QCs), the Monte Carlo approach is more suitable (natu- 
ral) for modeling the readout from a QC. 

III. ONE EXCESS ELECTRON IN THE CHAIN 

Let us consider first the simplest case of a chain of 
QDs initially doped with one excess electron. Since the 
Hamiltonian ^ preserves the number of electrons and 
their spin states, the system is restricted throughout its 
evolution to the one-electron Hilbert space Tii, and more 
precisely to its subspace pertaining to the initial spin 
state of the excess electron. Thus, the total wavefunction 
reads 

N 

\MT)) = Y.^n^)\j'^)' (3) 

where \ja) = a]a I*-*!' Ojv) denotes the state where one 
electron with spin a is located at the jth dot. The time- 
evolution of |'0i(r)) is governed by the Schrodinger equa- 
tion, which yields 

= e,A^ + + (4) 

where ioi = ^nn+i = and h — 1. Clearly, the two sets 
of these amplitude equations, pertaining to the spin up 
and spin down states, are completely equivalent and de- 
coupled from each other, i.e., the subspaces hX and 
are closed. As a result, if the excess electron is initially 
prepared in an arbitrary superposition of spin up and spin 

down states, |V'i(0)) = lit) + bi)' the two parts of 
the wavefunction will evolve symmetrically and indepen- 
dently of each other. This assertion is valid as long as 
one suppresses all uncontrollable spin-flip processes due 
to, e.g., presence of stray magnetic fields or spin-phonon 
coupling, which would otherwise destroy any pure spin 
state of the electron. Note that if the electrons in the ar- 
ray are to be used as qubits or carriers of quantum infor- 
mation, preserving their spin-coherence over long times 
becomes vital. This may potentially be realized by ap- 
plying a strong transverse magnetic field, that would re- 
move the spin degeneracy and make the spin-fiip process 
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FIG. 2: (a), (b) Single-electron transport in a chain of A'^ = 20 
QDs with no disorder and dissipation {Se = St — "/ = 0) for 
the case of (a) equal couplings, tjj+i = to, and (b) optimal 
couplings tjj+i = to^J {N — j)j, between the dots . (c),(d) 
Monte-Carlo simulations of detector response averaged over 
5000 trajectories for the chain with disorder and dissipation 
{Se = 0.1, St — 0.05 and 7 = 0.2) for the case of (a) equal 
couplings, and (b) optimal couplings between the dots. All 
parameters are normalized by to and the time is in units of 



energetically unfavorable. Experimental measurement of 
spin-relaxation times indicate sub-MHz rates |42| . 

In Fig. [2Ia) we plot the time-evolution of the occu- 
pation probabilities of TV = 20 QDs along a uniform 
(sj = SQ,tij = to), isolated from the environment chain, 
for the initial state |-!/'i(0)) — The initially well lo- 

calized electron wavepacket seen to propagate along the 
chain and spread to a larger number of dots, while si- 
multaneously splitting into several smaller components, 
which is due to the asymmetry of the initial conditions 
with respect to the electron forward and backward prop- 
agation. When the leading edge of the wavepacket is 
reflected from the end of the chain, its forward and back- 
ward propagating components begin to interfere with 
each other and after several reflections from the bound- 
aries of the chain, this interference causes the complete 
delocalization of the wavepacket over the entire chain. 

The results shown in Figs. [21 a), (b) have been obtained 
through the numerical solution of Eqs. These solu- 
tions, however, can be expressed in an analytical form ob- 
tainable through the Laplace transform. If one eliminates 



the first term via the transformation A° 



which is equivalent to the interaction picture, taking the 
Laplace transform of the set of Eqs. leads to 

isKf{\) - zA^"(0) = ^oi^^l(A) + ioif^i(A), 

where we have assumed equal energies £0 and tunnel- 
ing rates t^. Considering the case of only dot 1 oc- 
cupied initially, which means ^"(0) = 5ij, the quanti- 
ties if" are expressed in terms of the determinants of 
the coefficients. For N dots, the determinant I?jv(A) is 



readily shown to obey the recursion relation 2? at (A) = 
AI?Ar_i(A) + tlVN-2{X) with Vo = l and Vi ^ A, which 
identifies it with the Chebyshev polynomial of the second 
kind I?Ar(A) — n^i(A — A^) known to have the roots 

Xk — 2ito cos (y^ ^^i j ■ The determinants in the numera- 
tors of the expressions for _fC"(A) have similar properties. 
Taking advantage of this connection with the Chebyshev 
polynomials, one can, after some algebra, write the solu- 
tions for the amplitudes as 
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It is thus evident that the eigenstates of the coupled 
system oscillate with incommensurate frequencies corre- 
sponding to the roots Xk of Pjv, which in fact become in- 
creasingly densely spaced with increasing N. As a conse- 
quence, the system never revives fully to its initial state. 
This result has also been obtained in an entirely different 
context of driven multilevel systems H, 0, Q • 

Clearly, it is highly desirable to tailor the parameters 
of the system so as to achieve a non-dispersive trans- 
fer of the single-electron wavepacket between the two 
ends of the chain. This can be accomplished through 
the judicious choice of the interdot tunneling matrix el- 
ements, namely by choosing tjj+i according to tjj+i = 

to ^/{N-j)j, J = 1 , . . . , TV- 1 [Fig. Elb)] , to be referred to 
hereafter as optimal coupling. Then again, by exploring 
the properties of the Jacobi polynomials, one can obtain 
an analytic expression for the amplitudes, which in fact 
is even simpler, given by the binomial form 



A^- 1 



1/2 



-I sin (tor)] (^-1) cos (tor)^^~^\ 



while the eigenstates of the system have commensu- 
rate frequencies A^ = to(2fc — A^ — 1), and the electron 
wavepacket oscillates in a perfectly periodic way between 
the first and the last dots, whose occupation probabili- 
ties are given, respectively, by \A"\^ = cos (to'r)^*'^"^^ 

and |^5vP ~ si'^ (^o''")^^^"^''; a behavior reminiscent of a 
two-state system, in a generalized sense. 

Having established the conditions for the nondisper- 
sive transfer of the one-electron wavepacket between the 
ends of the chain, the next issue is the time-scale of the 
persistence of the coherence in the presence of disorder 
and/or fluctuations. In an actual experiment, these ef- 
fects will also be reflected in the signal detection, i.e., 
the measurement by a detector connected to one end of 
the chain. As already mentioned in the previous Section, 
this process is conveniently modeled through the quan- 
tum Monte Carlo technique. Although not necessary for 
the main issue here, we include it for the sake of com- 
pleteness, as it represents a mechanism of dissipation in- 
evitably present in an experiment. Continuing thus with 
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the case of the chain initially doped with one excess elec- 
tron, \tpi{0)) = I Iq), we plot in Figs. |2Ic),(d) the results 
of Monte-Carlo simulations for the detector signal aver- 
aged over 5000 trajectories. We chose the amplitudes 
of the energy-level and tunneling-rate fluctuations to be 
Se — O.lto and St = 0.05io, respectively, which, for typ- 
ical experimental parameters |4l| correspond to Se ^ 5 
jjLeV and St ^ 2.5 /xeV. The decay rate 7 of the electron 
from the last dot has been taken to be 7 = 0.2^0 ^ 2.4 
GHz. In Fig. |2Ic) we show the average detector signal 
for the case of equal mean couplings between the dots, 
while in Fig. El^d) the case of optimal mean couplings 
is illustrated. Comparing these figures with those corre- 
sponding to the isolated uniform chain [Figs. [21[a),(b)], 
as expected, the detector signal is found peaked around 
times when the electron occupies the last dot. Over the 
time scale of t '--^ Se~^,St~^ the system decoheres sig- 
nificantly due to the inhomogeneity introduced by the 
energy and coupling fluctuations, while the total occupa- 
tion probability of the chain decays roughly according to 
~ exp {—jt/N); the latter being the signature 
of dissipation. 

IV. TWO EXCESS ELECTRONS IN THE CHAIN 

We consider next a chain of QDs that is initially doped 
with two excess electrons. In analogy with the one- 
electron problem, the evolution of the system is restricted 
to the two-electron sector of the Hilbert space 7^2^'. We 
assume that the intradot Coulomb repulsion is so large, 
U ^ £j, tij , that the two electrons are forced to occupy 
different dots and thus the states Ijajjp) corresponding 
to two electrons at the same dot j are practically forbid- 
den. Then the state-vector of the system at an arbitrary 
time T is given by 

N 

where \ia,ja) = aL'^]/3 l^i' '^Jv) denotes the state 
where two electrons with spins a and (3 {a, f3 € {T,i}) 
are located at the ith and the jth dots, respectively. The 
corresponding amplitudes obey the following equations 
of motion 

= is. + e, + V,,)B'^+t,.uB:^,^+t,,+,B'^,^ 

and obviously, the four subspaces {hI^ Ti-l^ , 7^2^ and 
■^2^) of 7i2^ are governed by equivalent sets of equa- 
tions. These sets are decoupled from each other, since 
the Hamiltonian does not contain terms that directly 
couple different spin states. Note that since the states 
ba, j/3) are excluded from state-vector the transient 
Heisenberg coupling Jq = 4iQ/J7 between the spins is 
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FIG. 3: Two-electron transport in a uniform {Se — St — 0), 
isolated (7 — 0) chain of A'^ = 20 QDs for the initial state 
\tp2iO)) = |la,2;3). (a) Equal couplings between the dots, 
tjj+i = to, and (b) optimal couplings tjj+i = to\/{N — j)j 
{j = 1, A'^ — 1) and no interdot repulsion, V = 0. (c) 
Equal couplings and weak repulsion V = 0.75, and (d) equal 
couplings and strong repulsion V = 2.5. Note the different 
scales of the time axis. 



not accounted for in our model. This is justified pro- 
vided during the time r over which we consider the sys- 
tem dynamics the probability of this second-order spin- 
exchange process remains small, Jqt <^ 1, which in turn 
requires that U ^ 4tgr Then, due to the linearity of 
the Schrodinger equation, throughout the evolution each 
electron will preserve its spin state and the two electrons 
will not penetrate through one another, remaining thus 
in principle distinguishable. 

Unlike the one-electron case, useful analytic solutions 
do not seem easy to come by in the presence of a second 
electron, except for special cases. Also the analogy with 
the driven multilevel systems exists no longer. Obtained 
through the numerical solutions of Eqs. ©, in Fig. |2| 
we plot the time-evolution of the occupation probabili- 
ties of the QDs along a uniform, isolated chain for the 
initial state |?/'2(0)) = |lct)2^). Consider first the case of 
vanishing interdot Coulomb repulsion V — 0. Similarly 
to the single-electron case, the two-electron wavepacket 
propagates along the chain, spreads to a larger number 
of dots and, after several refiections from the boundaries, 
becomes delocalized over the entire chain, due to the 
multiple interference of its forward and backward propa- 
gating components [FiglS^a)]. These interference effects 
can be compensated for by employing the optimal cou- 
pling between the dots. Then every time the wavepacket 
reaches the end of the chain, the two electrons become 
fully localized at the last two dots [Fig|3Jb)] . The system 
is then formally (mathematically) equivalent to a chain 
of 2iV — 3 QDs doped with one electron. The solution for 
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FIG. 4: Energy-level diagram for two-electron states in the 
chain of QDs. 



the amplitudes has the following analytic form 



(j~if{N~l)\{N-2)\ 



1/2 



(r-l)\{j^l)\{N-i)\{N-j)\\ 
-I sin (toT)]^+^-3 cos (tor)2(^-2)-(^+^-3» . 



while the system has 2N — 3 distinct, commensurate 
eigenstates = to(2fc - 2iV + 2), where 1 < A; < 27V -3. 

Consider now the case of a finite near-neighbor in- 
terdot repulsion y > and equal tunneling rates be- 
tween the dots. The magnitude of this repulsion deter- 
mines the energy mismatch between the states \ia,if3) 
corresponding to two electrons occupying adjacent dots, 
j — i = 1, and those corresponding to two electrons 
separated by one or more empty dots, j — i > 1. As 
we show in Fig. |SJc), for relatively small interdot re- 
pulsion V = 0.75ioj the spreading and delocalization 
of the two-electron wavepacket is accelerated, which is 
due to inhomogeneity of the eigenenergies of states in- 
volved in the system's evolution. For larger interdot re- 
pulsion V = 2.5to, however, we find a rather counter- 
intuitive behavior of the system: the propagation dy- 
namics of the two-electron wavepacket and its delocal- 
ization is slowed down [Fig. Eld)]. This stems from 
the fact that when the two electrons are localized at ad- 
jacent dots, the corresponding states \ja, {j + l)/3) are 
coupled directly only to the states \{j — 1)q, (j + l)/3) 
and \ja, {j + 2)/3) [Fig. But if the interdot repulsion 
is larger than the tunneling rate, V > to, these transi- 
tions are nonresonant and the two electrons appear to 
be bonded with each other. Then, during the system 
evolution, the latter states are only virtually populated, 
while the second-order process responsible for the direct 
transition \ja, {j + l)/3) \{j + 1)q, (j + 2)^) is reso- 
nant. The effective coupling strength of this transition is 
given by = tg/V < to, and therefore the two-electron 
wavepacket propagation is slowed down. Since the en- 
ergy of the two-electron bonded state is larger than that 
of two separate electrons, the bonded state is not stable 
and non-adiabatic perturbations, such as collisions with 
the chain boundaries, gradually destroy it [Fig. E^d)]. 

Consider next the time-evolution of the system for 
the case when the two electrons are initially localized 
at the opposite ends of the chain, \tp2{0)) = jla,-^/?) 
(Fig. EJ. Upon propagation, the electrons collide with 
each other in the middle of the chain and then sepa- 
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FIG. 5: Two-electron transport in a uniform, isolated chain 
of iV = 20 QDs for the initial state |V'2(0)) = jl<:«,20,3>- (a) 
Equal couplings between the dots, tjj+i = to, and (b) op- 
timal couplings tjj+i = to \J (N — j)j {j — 1, N ~ 1) and 
no interdot repulsion, V = 0. (c) Equal couplings and weak 
repulsion V = 0.75, and (d) equal couplings and strong repul- 
sion V = 2.5. 



rate again. For each electron, the presence of the other 
represents an effective time-dependent potential barrier, 
which restricts the number of available dots per electron. 
Therefore, while after each collision the system partially 
revives to its initial state, the overall spreading and de- 
localization of the two wavepackets is slightly enhanced 
as compared to the case of Fig. E| Again, choosing the 
interdot tunneling rates according to the optimal rates 
tjj-i-i = ~ j)j , one can compensate for the inter- 

ference effects and achieve a perfect periodicity in the 
system dynamics [Fig. El^b)]. Finally, we note that in- 
creasing the interdot Coulomb repulsion strength results 
in acceleration of delocalization even further and inhibi- 
tion of occupation of the central dots of the chain where 
the two electrons collide [Fig. I3c),(d)]. 

We proceed now to the application of the quantum 
jumps method to the two-electron problem. Before the 
first jump occurs, the evolution of the system is restricted 
to the two-electron Hilbert space 712^. Immediately after 
the jump, the system is projected onto the single-electron 
Hilbert space Hf where it evolves until the second jump 
event. The total wavefunction of the system is thus given 
by 1^) = |'02) + jV'i) + jo), where IV'1,2) are defined in 
Eqs. (EJ and lO, respectively. In Fig. El we plot the 
results of our Monte Carlo simulations for the detector 
signal, averaged over 5000 trajectories. Due to limita- 
tions on the computer CPU time, we simulate a chain of 

= 10 QDs, in which the energy and tunneling-rate fluc- 
tuations and the decay rate from the last dot are taken to 
be the same as in Fig. EIc),(d), Se = O.lio, St = 0.05io 
and 7 = 0.2^0- (We do reaHze that TV = 10 - 20 may 
be demanding for present day technology, which, how- 
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ever, continues evolving.) In the two-electron case, the 
dynamics of the system is more involved as compared 
to the single-electron case, and the interpretation of the 
behavior of the detector signal is more subtle. At small 
times, while the probability that the first jump has oc- 
curred is not yet very large ||V'2(t)|P ^ 1, one can vi- 
sualize the behavior of the detector signal by compar- 
ing the four plots in Fig. with the corresponding 3D 
graphs in Figs. |3|and|Sl One thus sees that the signal 
is peaked around times when the population of the last 
dot is maximal. At larger times, however, the averaged 
over many trajectories signal represents an interplay be- 
tween the two-electron propagation before the first jump, 
and single-electron propagation after the jump. Since the 
jumps occur at random times, the single-electron prob- 
lem at each trajectory has different initial conditions, and 
therefore at larger times the detector signal does not have 
any pronounced structure except for an overall exponen- 
tial decay with the rate 7/iV. Thus in Fig. El^a) the 
three pronounced peaks of the signal around r ~ 7, 18 
and 30 xt^^ correspond to the arrival of the two-electron 
wavepacket at the last dot. In Fig. Elb), where we simu- 
late the case of two-electron bonding via strong interdot 
Coulomb repulsion V = 2.5tQ, the first small peak of 
the signal around r ~ 7 x ^ corresponds to the small 
unbonded fraction of the two-electron wavepacket, while 
the second larger peak represents the delayed arrival of 

the bonded pair of electrons at the last dot (recall that 

(2) 

the effective tunneling rate t^g for the bonded electron 
pair is smaller than the single-electron tunneling rate to). 
Our numerical simulations show that increasing the in- 
terdot repulsion does indeed strengthen the bonding and 
reduces the height of the first peak. Next, in Fig. El^c), 
where we simulate the case of optimal couplings between 
the dots, the periodic behavior of the signal at short times 
reflects the periodic arrival of two-electron wavepacket at 
the last dot. Finally, in Fig. Eld), we illustrate the results 
for the case when the two electrons are initially localized 
at the opposite ends of the chain. Here, before the first 
jump has occurred, the second electron is effectively con- 
fined in a chain of length N/2 QDs. As a result, the signal 
initially oscillates with twice the frequency as compared 
to the case of Fig. EJa), while after the jump, the entire 
chain becomes accessible to the single remaining electron. 



V. ENTANGLEMENT OF TWO-ELECTRONS 

Material particles such as electrons represent a promis- 
ing alternative to circumvent communication and detec- 
tion loopholes that occur during the tests of Bell inequal- 
ities l4q . Despite an ever growing number of theoretical 
proposals 46, 47, 48, 49], the experimental demonstra- 
tion of entanglement generation and detection in con- 
densed matter systems is still an open issue. One of the 
hardest steps towards such a realization is the necessity 
of controllable spatial separation of the two constituents 
of an entangled pair [48| . A similar problem would hin- 
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FIG. 6: Monte-Carlo simulations of detector response for the 
two-electron transport in a chain of A'^ = 10 QDs with dis- 
order and dissipation {5e = 0.1, St = 0.05 and 7 = 0.2), 
averaged over 5000 trajectories. In (a), (b) and (c) the ini- 
tial state is ji/'2(0)) — |1q,2^), and in (d) the initial state 
is |7/'2(0)) — |la,10/3). Other parameters are: (a), (d) equal 
mean couplings between the dots and no interdot repulsion, 
V = 0; (b) equal couplings and strong repulsion, V — 2.5; (c) 
optimal couplings and no repulsion, I^ = 0. 



der the realization of large-scale solid-state QCs. One of 
the main difficulties with the existing proposals for inte- 
grated QD based QCs 25, 26, 29, 30J is that the qubits 
(electron spins) interact with the nearest neighbors only, 
and there is no efficient way of transferring the infor- 
mation between distant qubits (one can swap the qubit 
state with its neighbor qubit, then the latter with the 
next neighbor, etc., which is inefficient since each SWAP 
action requires precise dynamical control of the parame- 
ters) . Our intention in this section is to show that use of 
the optimal coupling circumvents these difficulties. 

We adopt here an entanglement scheme (entangler) 
proposed in [s^l , which is based on the exchange in- 
teraction in a double-dot system where each QD is oc- 
cupied by a single electron. When the interdot potential 
barrier is high, the tunneling is negligible, te ~ 0, and the 
electronic orbitals are well localized at the corresponding 
dots. By lowering the interdot tunnel barrier, one can 
induce a finite overlap between the orbitals. Then the 
two electrons of this "artificial molecule" will be subject 
to a transient Heisenberg coupling [2^ [s^, E3 



H^ir) ^ ^JiT)SL ■ Si, 



(7) 



where Sl,r. are the spin— 1/2 operators for the left (L) 
and right (i?) QDs, respectively, and J(r) — At1{T)/U 
is the effective time-dependent spin exchange constant. 
If initially the two electrons have opposite spin states, 
= \L^,Ri), and the pulsed Heisenberg couphng is 
applied for a specific duration such that = J J{T)dT = 
7r/2, the square root of swap action (-y/swAp) is realized. 
Then the resulting state of the system is the maximally 
entangled state |(/)) — > ^( i?|) -I- i |L|, i?|)). 

Thus, a possible scheme, capable of producing meso- 
scopically separated EPR correlated pairs of electrons. 
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(a) 



Entangler 




FIG. 7: Three step generation of mesoscopically separated 
entangled electron pair, (a) In step 1, employing the optimal 
coupling for the chains 1 — L and R~ N, the two electrons are 
transferred to the entangler. Raising the barriers Br and Bl, 
the entangler is temporarily isolated and the two electrons are 
trapped in dots L and R, respectively, (b) In step 2 a pulsed 
Heisenberg exchange coupling is applied, by turning on and 
off the tunneling matrix element t^{T). (c) In step 3 the two 
electrons are transferred back to dots 1 and A''. 




FIG. 8: Monte Carlo simulation of entangled state gener- 
ation in a chain of A'^ = 20 QDs. The initial state is 
\i'2{Q)) ~ I'^u^l)- The occupation probabilities of the 
two-electron states ll|,A'^|) (thin solid line), |L|,J?|) (thin 
dashed line), -^( , + i\Li, Rf)) (thick dashed line), 
and , Ni) + i\li,N-j)) (thick solid line) are plotted 

as functions of time for the three-step process illustrated in 
Figs. IZIa)-(c), respectively. The parameters are f™'"' = 6to, 
U — lOOto (Ar ~ 0.551^^) and the optimal coupling is em- 
ployed separately for each subchain 1 — L and R — A^. Other 
parameters are the same as in Fig. |H| 



relies on the double-dot entangler that is embedded in a 
long chain of TV = 2M coupled QDs. To be specific, let 
us consider a symmetric structure depicted in Fig. CJa), 
in which the entangler, defined by the barriers and 
Bfj, is represented by the two middle dots of the chain, 
L = M and R — M + 1. As discussed at the end of 
this Section, the exact location of the entangler is not re- 
strictive for the success of the scheme. The entanglement 
scenario then proceeds in three steps illustrated in Figs. 
[7fa)-(c), respectively. During the first step [Fig. [7{a)], 
the potential barrier between the dots L and R is kept 
high and the exchange coupling is closed, = 0. There- 
fore the two parts of the chain, composed of dots 1 — M 
and (M -I- 1) — N, respectively, are decoupled from each 
other. The two unentangled electrons, initially located 
at dots 1 and N, are coherently transfered to the en- 
tangler by employing the optimal coupling in both parts 
of the chain. At time r — 7r/(2io), when each electron 
reaches the corresponding dot of the entangler, the bar- 
riers Bl and Bj^ are pulsed to higher potentials so as to 
trap the corresponding electron. Note that, since = 0, 
the trapping process need not occur simultaneously for 
both electrons. At the second step [Fig. [3b)], the ex- 
change interaction is applied, by adiabatically turning 
on and off the tunneling matrix element according to 
t^(r) = tf^^'sechlir - T'"^'^)/Ar], where C'^'^ is the peak 
amplitude of the pulse and Ar is its width. These param- 
eters must be chosen such that the pulse area 9 satisfies 
9 = n/2, which yields {t'^^^' fAT = nU/ld. At the same 
time, according to the adiabaticity criteria '2^, the con- 
ditions AT~^,t™'^'^ <C Ae, U should be satisfied. Finally, 
at the third step [Fig. EJc)], the barriers and Bj^ are 
reset to their initial values and the two electrons propa- 



gate along the chain and after time r = 7r/(2to) become 
fully localized at the opposite ends of the chain, at dots 
1 and N. 

If the two electrons are prepared initially in the oppo- 
site spin states, [^'2) = |lT;^i)i during the three-step 
process the system will evolve into the maximally en- 
tangled state IV2) ^ ^(|1|,7V4) +i\li,N^)). If, on 
the other hand, the two electrons are prepared initially 
in the same spin state, |l|,iV|) or |l|,iV|), this state 
will remain unchanged (to within the phase factor e'^^^). 
Finally we note that if, during the second step, the ex- 
change pulse area 9 is equal to tt (instead of 7r/2), the two 
electrons will swap their states |1q, Np) — > z Na). 

We have performed Monte Carlo simulations of the 
above scenario for an even total number of nearly iden- 
tical QDs and the initial state 1-02) = |1|,A''|). The 
evolution of the occupation probabilities \{<pi\4'2{T))\'^ 
for the two-electron states \(j)o) — |l|,Af|), \<j)i) — 
\L^,Rl), 102) = -^{m,Ri)+i\L^,R^)), and [(/.g) = 
-^( |l|,iV|) + i |l|,iV|)) during the three stages of the 
process are plotted in Fig. |S1 One can see that the fi- 
delity of the entangled state generation at the end of 
the process is very high, |(03|'02(T'cnd))|^ — 0.98. Em- 
ploying spin-sensitive single-electron transistors [i^, ^3 ' 
one could then use this EPR pair to perform a test of 
Bell inequalities with material particles (electrons). The 
controlled two-electron coherent dynamics studied above 
could be also envisioned in an integrated quantum regis- 
ter. Such a register would be composed of a large num- 
ber of sub-registers, each containing two or more adja- 
cent qubits represented by spins of single electrons in 
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individual QDs. The sub-registers are embedded in a 
two-dimensional array of empty QDs. In combination 
with the mechanism of controlling the tunnel-coupling 
between the dots, this two-dimensional grid would real- 
ize a flexible quantum channel, capable of connecting any 
pair of qubits within the register. Thus, to transfer the 
information, one connects distant sub-registers by a chain 
of empty QDs and arranges the optimal tunnel-coupling 
between the dots to achieve a non-dispersive propaga- 
tion of the qubit, followed by its controlled entanglement 
or SWAP with a target qubit. Note that this scheme is 
analogous to a proposal for an integrated ion trap based 
QC |5^]. There, in order to circumvent the difficulties 
associated with single large ion trap quantum register, it 
was proposed to use many small sub-registers, each con- 
taining only a few ions, and connect these sub-registers 
to each other via controlled qubit (ion) transfer to the 
interaction region (entangler) represented by yet another 
ion trap. It would seem that contemplating subregisters 
consisting of few (as few as two) QDs is no less practical. 

Before closing this section, let us address the issue of 
location of the entangler. Since, as was noted above, the 
two electrons need not arrive to the entangler simultane- 
ously, it is not necessary to position the entangler in the 
middle of the chain or synchronize the propagation of the 
two electrons during the first and the last steps. Actu- 
ally, the entangler can be realized by any pair of dots in 
the chain, and in particular, even by the first two dots. 
Then only the second electron, initially located at dot N 
should be brought into the entangler to interact with the 
first electron. 



VI. CONCLUSIONS 

To conclude, we have studied the microscopic dynam- 
ics of one or two electron transport in a linear array of 



tunnel-coupled quantum dots. We have shown that co- 
herent electron wavepacket propagation in this system 
reveals a number of novel, surprising effects, such as 
suppression or enhancement of interference and spread- 
ing, two-electron bonding and collision. We should note 
that the coherent electron wavepacket propagation in ar- 
rays of tunnel- coupled QDs bears many analogies with 
spin- wave propagation in spin chains |51| . electromag- 
netic wave propagation in nonlinear periodic media |52| 
and matter-wave propagation in optical lattices, where 
a transition from superfluid to Mott-insulator phase has 
recently been observed 0, H- With an unprecedented 
control over system parameters, arrays of QDs allow for 
studies of numerous coherence and correlation effects in 
many-body physics. As an example, dynamical localiza- 
tion of electrons is expected to be observed in the ar- 
rays with large disorder fc, 5t ^ to 0,^3- Our studies 
are also aimed at possible applications in future quan- 
tum computation schemes (23, i2S | , as many solid state 
quantum computer pro pos als involve a rray s of interact- 
ing quantum dots [lllllllllMlllllIllllll. We 
stress that if the electrons in the array are to be used 
as qubits or carriers of quantum information, preserv- 
ing their spin-coherence over long times becomes vital. 
Two-dimensional arrays of electrostatically- and tunnel- 
interacting quantum dots have also been suggested as 
an efficient tool for analog computation (Sj- In addi- 
tion, there is presently growing interest in implementing a 
quantum cellular automata chain '53| , in which the com- 
putation is performed by applying a sequence of unitary 
transformations to the input state, upon its propagation 
through the chain towards the output. Such a chain with 
tunable couplings between its links would correspond to 
a universal quantum computer, closely resembling the pi- 
oneering idea of Feynman [55| to simulate one physical 
system by another. 
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